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In this work, we have studied accretion of the dark energies like new variable modified Chaplygin 
gas (NVMCG) and generalized cosmic Chaplygin gas (GCCG) onto Schwarzschild and Kerr-Newman 
Black holes. We find the expression of the critical four velocity component which gradually decreases 
for the fluid flow towards the Schwarzschild as well as Kerr-Newman Black hole. We also find 
the expression for change of masses of the black hole in both cases. For the Kerr-Newman black 
hole which is rotating and charged we calculate the specific angular momentum and total angular 
momentum. We showed that in both cases due to accretion of the dark energy mass of the black 
hole increases and angular momentum increases in case of Kerr-Newman black hole. 
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I. INTRODUCTION 



Different observational data together with observations of supernovae of type la [H-Q) WMAP 
Chandra X-ray Observatory [6| strongly indicate that our universe is undergoing an accelerating phase. 
Nonbaryonic matter recognized as dark energy having negative pressure and violate the strong energy 
condition may explain this accelerated expansion. There are various candidates to play the role of the dark 
energy which is the dominant part of the universe. The dark energy candidates are cosmological constant 
[ioj . quintessence, K-essenc e fill . Chaplygin gas [l2|, its modification known as modified Chaplygin gas 
(MCG) QJ], tachyonic field QjjTDBI-essence etc. The Chaplygin Gas (EoS p = — — , B > 0) 12] acts 
as pressureless fluid for small value of the scale factor and tends to accelerated expansion for large value 
of scale factor. Generalization of Chaplygin gas model known as Generalized Chaplygin gas which satisfies 
p=--^,0<a<l tl6Hl8j- This model also modified to Modified Chaplygin Gas (MCG) having the 
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Eos v = Ap - -pL < a < 1. A > 0. B > [13l Il9l |20|. That illustrates a radiation era (A = 1/3) while 
the scale factor is vanishingly small and ACDM model for infinitely large scale factor. Further Guo and 



Zhang 21 1 established Variable Chaplygin Gas with Eos is p = — — , where B = B(R), R is the scale factor 
and B is a positive function of the scale factor. Subsequently Debnath [221 ] provided Variable modified 
Chaplygin Gas with Eos is p = Ap — for the accelerating phase of the universe. The another candidate 
of dark energy was introduced by Chakraborty et al (23j . known as New Variable modified Chaplygin 
Gas (NVMCG) which follows the equation p = A(R) — B ^ , < a < 1 which gives interesting physical 
significance. In 2003 P.F. Gonzdlez-Diaz [24| gives the idea of another form of dark energy to the conse- 
quence of accelerating phase of universe namely Generalized Cosmic Chaplygin Gas (GCCG), This model is 
stable and free from unphysical behaviour even when the vacuum fluid satisfies the phantom energy condition. 

A cosmolo gica l property in which there is an infinite expansion in scale factor in a finite time termed 
as 'Big Rip' [26|. In the phantom cosmology, big rip is a kind of future singularity in which the energy 
density of phantom energy (p + p < 0) will become infinite in a finite time. To realize the Big Rip scenario 
the condition p + p < alone is not sufficient [27[ • Distinct data on supernovas showed that the presence 
of phantom energy with —1.2 < w < —1 in the Universe is highly likely [28j |. The accretion of phantom 
dark energy onto a Schwarzschild black hole was first modelled by Babichev et al [29] . They established 
that black hole mass will gradually decrease due to strong negative pressure of phantom energy and finally 
all the masses tend to zero near the big rip where it will disappear. Accretion of phantom like modified 
variable Chaplygin gas onto Schwarzschild black hole was studied by Jamil [3(| who showed that mass of 
the black hole will decrease when accreting fluid violates the dominant energy condition and otherwise will 
increase. Also the accretion of dark energy with EoS p — wp onto the Kerr-Newman black hole was studied 
by Madrid et al [3l| and they obtained that if w > —1, mass and angular momentum increase. Mass of the 
black hole grows up unboundedly whereas the angular momentum increases up to a given level. 

In the present work, we have studied accretion of dark energy namely new variable modified Chaplygin gas 
(NVMCG) and generalized cosmic Chaplygin gas (GCCG) onto Schwarzschild as well as most generalized 
Kerr-Newman black holes. For natures of black hole mass function with angular momentum have been 
analyzed when NVMCG and GCCG like dark energies accrete upon black holes. 



II. ACCRETION OF DARK ENERGY ONTO SCHWARZSCHILD BLACK HOLE 

Let us consider a spherically symmetrical accretion of the dark energy onto the black hole. We consider a 
Schwarzschild black hole (static) of mass M which is gravitationally isolated (in geometrical units, G = 1 = c) 
[2H [30( 1 described by the line element 



ds 2 



(l - ^ dt 2 -(l- 1 dr 2 - r 2 (d9 2 + sin 2 Odcf) (1) 



where, r being the radial coordinate. Energy momentum-tensor for the DE, considering in the form of 
perfect fluid having the EoS p — p(p), is 

= (p + p)UftU v - pg^u (2) 

where p, p are the density and pressure of the dark energy respectively and u M = is the fluid 4- velocity 
satisfying u^u^ = 1. We assume that the in- falling dark energy fluid does not disturb the spherical 
symmetry of the black hole. 

The relativistic Bernoulli's equation after the time component of the energy-momentum conservation law 
T.^f — provide the first integral of motion for stationary, spherically symmetric accretion onto BH which 
yields 

(p + p) (l- ^+u 2 Y x 2 u = d (3) 
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where x = -j? and u = Of- is the radial component of the velocity four vector and C\ being the integrating 
constant. In the case of fluid flow directed towards the black hole, we must have u < 0. 

Moreover, the second integration of motion is obtained from u^T.^f = 0, which gives 



ux exp 



dp 



P + p(p) 



-A 



(4) 



where, A(> 0) is the integration constant, poo is the dark energy density at infinity. Further value of the 
constant A is evaluated for different DE model. 



Using the equations (3) and (4) we get 



(p + p) \ 1 Yu 2 exp 



dp 



Poo P + P{p) 



C 2 



(5) 



where C2 = —C\/A = poo + p(poo)- 

If n be the concentration of dark energy which satisfies the following equations 



n{p) 



exp 



dp 



> Pao P + p{p) 

where rioo = n(/9oo) being the concentration of the dark energy at infinity. 

The constant value A can be determined by finding the critical point of the accretion using [32| then, 



(G) 



I 



C 2 s(P*) 



l-3u 2 



(7) 



where c s = \/ ^ is the usual speed of sound, it, is the critical four velocity component and is the 



density at the critical point x*. One may noted that for c 2 s > or c 2 s < 1, no critical point exists outside the 
black hole (i.e., x* > 2). Using (5) and (7), we get the following relation 



p* +p(p*) 

Poo +p(Poo) 



= yl + 3c 2 (p*) exp 



dp 



P + p(p) 



(8) 



The rate of change of mass M of the black hole is computed by integrating the flux of the dark 
energy over the entire horizon of the black hole i.e., M = § T[dS, where T[ represents the radial compo- 
nent of the energy momentum densities and the surface element of the black hole horizon dS = ^ —gdQdij) ^29] . 



Using the above equations we obtain the rate of change of mass as 



M = An AM 2 (p + p) (9) 

Since the Schwarszchild black hole is static, so the mass of the black hole depends on r only. When some 
fluid accretes outside Schwarszchild the black hole, the mass function M of the black hole is considered as 
a dynamical mass function and hence it should be a function of time also. So M of the equation (9) is time 
dependent and the increasing or decreasing of the black hole mass M sensitively depends on the nature of 
the fluid which accretes upon the black hole. 

At the black hole horizon (r = 2M i.e., x = 2) the relation between four velocity ur = u(pn) and the 
energy density pu at the black hole event horizon is given by 
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A p H +p{p H ) 
4 Poo +p{poo\ 



A 2 



16u% 



exp 



dp 



P + p{p) 



and velocity four component at the horizon of black hole is given by 



u H 



■ exp 



dp 



p + p(p) 



(10) 



(11) 



A. Model I: New Variable Modified Chaplygin Gas as dark energy model 

We consider the background spacetime is spatially flat represented by the homogeneous and isotropic FRW 
model of the universe which is given by 



ds 2 = -dt 2 + R 2 (t) [dr 2 + r 2 {d9 



2 i -2 

1 sin 



'1] 



where R(t) is the scale factor. We assume the universe is filled with New Variable Modified Chaplygin 
Gas (NVMCG) and the EoS is [1| given by 



p = A'{R)p- 



B{R) 



with < a < 1 



(12) 



where A'(R), and B(R) are function of the scale factor R. For a particular choice A'(R) = AoR~ n and 
B(R) — BoR~ m with Aq, Bq, to, n are positive constants. For n = m = 0, this model reduces to modified 
Chaplygin Gas, and for n = the model reduces to the variable modified Chaplygin gas model. 

The Einstein's equations for FRW universe are (choosing G = c = 1) 



H 2 



8tt 



(13) 



H = - 8 i(p + p) 

Conservation equation satisfied by the dark energy model NVMCG is 



(14) 



p + 3H(p + p) =0 

where H = is the Hubble parameter. 

Expression for the energy density for NVMCG model is obtained from (12) and (15) as [2J 



(15) 



p = R 3 exp 



/3A R~ n 



Co + ^ X ^^ T ( m -^ 1 + a \ XR -n 

Aq V n 



V n 

where Co is an integration constant, T(s,t) is the upper incomplete gamma function and X = 

Following [32j , the critical values for this model are as follows 

B 



(16) 



A'pl 



A' + a- 
P 

+ 1 +aB 



a+l 



p* L (l+3A')+3aB 
p? +1 (l+3A')+3aB 
2 [A'p^ 1 + aB] 



(17) 
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FIG. 1: Accreted fluid velocity with respect to radial coordinate x for NVMCG onto Schwarzschild black hole. 
FIG. 2: Accreted fluid velocity with respect to radial coordinate x for GCCG onto Schwarzschild black hole. 

Also from equation (6) and (4), we obtain the ratio of number density of NVMCG near the horizon and 
at the infinity as in the following 



n{p) 



B'-p° 



B'-p% 



where p, 



(l+a)(l+A') ' 



B' 



l+A 



7 and 



1 


>? +1 (l + 3A') + 3aB" 


3/2 


B' 
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A'p« +1 + aB 




.B' 





Dark energy density pn at the event horizon of the black hole satisfies the following equation 



\Ph +1 


-B'~ 


2/i- 




^ Poo\ a 


U +1 


- B'_ 






KPh) 



(18) 



(19) 
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Also the energy density can be expressed in terms of x as 



1-- + — 

x I x 



A 2 ( B' 



„a+l 1 M 



B' 



2a r 



B' 



B' 



2(1-A») 



1 



(21) 



Accreted fluid velocity with respect to radial coordinate x for NVMCG onto Schwarzschild black 
hole are drawn in figure 1. The fluid velocity decreases as x increases, i.e., accreted fluid velocity 
is high near the black hole and velocity is low when the fluid is far from the black hole. The black 
hole mass increases for the NVMCG dark energy type accreted fluid as time goes on, which is shown in 
figure 3. Also relative density p/poo of accreted fluid increases when x increases from the black hole (figure 5). 



B. Model II: Generalized Cosmic Chaplygin Gas as the dark energy Model 

A new version of Chaplygin gas which is known as Generalized Cosmic Chaplygin Gas (GCCG) [ID, [2£ 
obeys the equation of state 



C+ (p 1+a -C) 



(22) 



where C 



A" 



1, A" takes either positive or negative constant, — I < w < and / > 1. The EOS 



reduces to that of current Chaplygin unified models for dark matter and dark energy in the limit w — > 



6 



and satisfies the conditions: (i) it becomes a de Sitter fluid at late time and when w = — 1, (ii) it reduces 
to p = wp in the limit that the Chaplygin parameter A" — > 0, (iii) it also reduces to the EOS of current 
Chaplygin unified dark matter models at high energy density and (iv) the evolution of density perturbations 
derived from the chosen EOS becomes free from the pathological behaviour of the matter power spectrum 
for physically reasonable values of the involved parameters at late time. This EOS shows dust era in the 
past and ACDM in the future. 

From the conservation equation we have the expression for energy density in the form [25[ 



P = 



C 



1 



B 



Jj3(l+a)(l+tu) 

Following [32j , the critical values are obtained in the form 



(23) 



where y = 



= p* 1_ 


- a y 


y 




P l +a + 


3y 


pl +a + 







Ca+(pl+ a -C) 1 w {pl +a (a + w + aw)-Ca) 



(24) 



From equations (6) and (4), we obtain the ratio of number density of GCCG near the horizon and at the 
infinity as in the following 



n(p) 



MpW) 

fl(P<x)h{Poa) 



where, 

fl(p)= [p~ a {C + (-C + p 1+a )- w } - p] , h{p)=p^ 1+w) {Cp-<*-p) w and 
Energy density can be expressed in terms of x as 



1 



(l + a)(l + w) 



(25) 



(26) 



where A = — 



_ _1 / pl + "+3y\ 3/2 n(p,) 



Ml 



x 4 n(p) 



\ flip) ' 


2(«/-l) 


[ hip) ' 




./l(Poo). 




.hipoo). 





(27) 



Accreted fluid velocity with respect to radial coordinate x for GCCG onto Schwarzschild black hole are 
drawn in figure 2. The fluid velocity decreases as x increases, i.e., accreted fluid velocity is high near the 
black hole and velocity is low when the fluid is far from the black hole. The black hole mass increases for 
the GCCG dark energy type accreted fluid as time goes on, which is shown in figure 4. Also relative density 
pi Poo of accreted fluid increases when x increases from the black hole (figure 6). 



III. GENERAL ACCRETION OF DARK ENERGY ONTO KERR-NEWMAN BLACK HOLE 



For the rotating, charged, stationary, axisymmetric black hole, let us consider Kerr-Newman space-time 
metric (considering G = c = 1) [3lj prescribed by the line element 
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FIG.3 FIG.4 
FIG. 3: Changes of the mass with respect to time of NVMCG onto Schwarzschild black hole. 
FIG. 4: Changes of the mass with respect to time of GCCG onto Schwarzschild black hole. 
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FIG. 5: Relative density of accreted fluid with respect to radial coordinate x for NVMCG onto Schwarzschild black 
hole. 



FIG. 6: Relative density of accreted fluid with respect to radial coordinate x for GCCG onto Schwarzschild black 
hole. 



ds 2 = 1 



Q 2 - 2Mr 



dt 2 



2a(2Mr - Q 2 ) sin 2 9 



r 2 + a 2 cos 2 ( 



I 2 i 2 2 

[r + a cos 



de 2 



dtd( ^ r 2 +a 2 + Q 2 - 2Mr 
(2Mr-Q 2 ) a 2 sin 2 9 
r 2 + a 2 cos 2 9 



dr 2 



(28) 
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where M and Q are respectively mass, electric charges of the black hole. Also a = J/M is the specific 
angular momentum per unit mass and J is the total angular momentum of the black hole. Since Kcrr- 
Newman metric is static, the time evolution induced by accretion will be taken into account by the time 
dependence of the scale factor entering the integrated conservation laws and the rate equations for mass and 
angular momentum. Also the energy tensor T^ v for dark energy satisfies the relation (2). The first integral 
of motion comes from the time component of the energy-momentum conservation law T£ v = which can be 
put in the following form for Kerr-Newman metric: 



d 
d9 



d_ 

dr 

{p + p) I 1 + 



(P + P) 1 + 



Q 2 - 2Mr \ dt dr 



ds ds 



2r 



a 2 cos 2 9 



(p + p) 1 + 



Q 2 - 2Mr \ dt dr 



r 2 + a 2 cos- 



\ dtde ' 




) ds ds 


+ 



cos6 2a 2 sin 9 cos 9\ . , / Q 2 - 2Mr 
(P + p) 1 + 



sin 9 



a 2 cos 2 9 



r z + a 2 cos 2 



2 ^ac2 n l ds 
= 



(29) 



In general, p and p are functions of t, r and 9, but in our considered dark energy models p and p depend 
only on time t. Keeping 9 constant and radial four velocity component «=^we get the expression for rate 
of change of mass due to accretion of dark energy [3l| 



M = 



Jn 



AnA M M 3 r ( J \ 

dS = arctan I — j [p + p) 



(30) 



r sin 



with r and J are constants. Which gives the expression 
JdM 



hi — 



M 



M M 3 r arctan (^) 



= 4nA M f (p + p)dt 



with 



A M = -j-p ( r2 + a 2 cos2 9 ) ex P 



dp 



(31) 



(32) 



> Poo P + P(p) 

When a fluid flow directed toward the black hole we have u < and Am > is dimensionless constant. 

Again keeping r as constant and we get the rate of change of angular momentum of Kerr-Newman black 
hole as [3l| 



a = - / rTUS = 



where dS = r 2 sin 



2TT 2 A a ar 2 (p + p) 
V r 2 + a 2 

, 9 constant. Equation (33) simplifies to 



(33) 



with 



A„ 



/ 5 da = 2ir 2 A a / (p + p)dt 

Ja n ar 2 J tn 



1 



w sin 9(r 2 + a 2 cos 2 9) exp 



dp 



p + p{p) 



constant. 



(34) 



(35) 



When a fluid flow is directed toward the black hole then u = W- < for A a > 0. 

Expression for the rate of change of total angular momentum, with M, r constant, is given by [3l| 



(MrT® + aTS) dS = tt( P + p) 



2JnA a r i A , , 2 (J 
• 4A M M 2 r arctan I 



\Mr 



(36) 
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Which reduces to 



-f 

J. J a 



dJ 



2JirA a r 



4^ M M 2 rarctan( 1 ^:) 



7T f ( 

Jt 



(p + p)dt 



We shall obtain the quantity M, a, J for the following cosmological models: 



(37) 



A. Model I: New Variable Modified Chaplygin Gas as dark energy model 

For the NVMCG satisfying the equations (11-15) we derive the following quantities 



8ttA 2 m 



hi = VPo - \/~P 



(38) 



2n 3 A 2 „ 



(39) 



i j = Vpo - Vp 



(40) 



where p satisfying the relation (16) and from this equation we obtain the present value of density p which 
is given by 



po = R exp 



3 , 3A R ' 



B f 3A (l + a) 
A \ n 



3(1 + q) + ^ 



f j m-3(l + a) 3^0(1 + 0;) ^. 



(41) 



where Ro is the present value of scale factor. Changes of the mass with respect to time of NVMCG onto 
Kerr-Newman black hole with the constant J and variable J are drawn in figures 7 and 8 and they are 
increasing. The angular momentum and the total angular momentum with time are shown in figures 9 and 
10 and they are increasing with time. 



B. Model II: Generalized Cosmic Chaplygin Gas as the dark energy Model 

For the Generalized Cosmic Chaplygin Gas (GCCG) satisfying the equations (12-14) and (21,22) we derive 



R 3(l+ a )(l+w) = 
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2(l+a) 
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2(l+a) 
filj) c 



(l+w) 



(42) 



(43) 



(44) 
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FIG. 7: Changes of the mass with respect to time of NVMCG onto Kerr-Newman black hole with the constant J. 



FIG. 8: Changes of the mass with respect to time of NVMCG onto Kerr-Newman black hole with J = J(t). 
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FIG. 9: Changes of the specific angular momentum (a) with respect to time of NVMCG onto Kerr-Newman black 
hole with the constant r. 

FIG. 10: Changes of the total angular momentum (J) with respect to time of NVMCG onto Kerr-Newman black 
hole. 



with, 



Po 



C 



B 

3(l+a)(l+ui) 




(45) 



11 




FIG.ll FIG.12 
FIG. 11: Changes of the mass with respect to time of GCCG onto Kerr-Newman black hole with the constant J. 



FIG. 12: Changes of the mass with respect to time of GCCG onto Kerr-Newman black hole with J = J(t). 




FIG.13 FIG.14 

FIG. 13: Changes of the specific angular momentum (a) with respect to time of GCCG onto Kerr-Newman black 
hole with the constant r. 

FIG. 14: Changes of the total angular momentum (J) with respect to time of GCCG onto Kerr-Newman black hole. 

Changes of the mass with respect to time of GCCG onto Kerr-Newman black hole with the constant J 
and variable J are drawn in figures 11 and 12 and they are increasing. The angular momentum and the 
total angular momentum with time are shown in figures 13 and 14 and they are increasing with time. 
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IV. DISCUSSIONS 

In this work, we have studied recently proposed two types of dark energy models like new variable 
modified Chaplygin gas (NVMCG) and generalized cosmic Chaplygin gas (GCCG). Accretion of NVMCG 
and GCCG onto the Schwarzschild and Kerr-Newman Black holes have been discussed. Our dark energy 
fluids violate the strong energy condition (p + 3p < in late epoch), but do not violate the weak energy 
condition (p + p > 0). So the models drive only quintessence scenario in late epoch, but do not generate 
the phantom epoch (in our choice). We find the expression of the critical four velocity component which 
gradually decreases for the fluid flow towards the Schwarzschild as well as Kerr-Newman Black holes. 
Astrophysically, mass of the black hole is a dynamical quantity, so the nature of the mass function is 
important in our black hole models for different dark energy filled universe. Previously Babichev et al 
[29| have shown that the mass of black hole decreases due to phantom energy accretion. We here found 
the expression for change of masses of the Schwarzschild and Kerr-Newman black holes in both the dark 
energy models and have seen that they are increasing in course of time. Since our considered dark energy 
candidates do not violate weak energy condition, so the dynamical mass of the black holes could not 
decaying by the accretion of dark energies, though the pressures of the dark energies are outside the black 
holes. The relative density p/poo increases as r increases outside the black hole. For the most generalized 
Kerr-Newman black hole (which is rotating and charged) we have obtained the specific angular momentum 
(a) and total angular momentum (J). We showed that in both cases due to accretion of the dark energy 
mass of the black hole increases and angular momentum increases in case of Kerr-Newman black hole. The 
mass of the black hole increases for constant and variable J. The relative density shows the same nature as 
well as Schwarzschild black hole. 
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